On the stability of Bose-Fermi mixtures 
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We consider the stability of a mixture of degenerate Bose and Fermi gases. Even though the bosons 
effectively repel each other the mixture can still collapse provided the Bose and Fermi gases attract 
each other strongly enough. For a given number of atoms and the strengths of the interactions 
between them we find the geometry of a maximally compact trap that supports the stable mixture. 
We compare a simple analytical estimation for the critical axial frequency of the trap with results 
based on the numerical solution of hydrodynamic equations for Bose-Fermi mixture. 
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Experimental realizations of mixtures of degenerate 
atomic gases have opened exceptional possibilities to in- 
vestigate fundamental many-body quantum phenomena. 
The reason is that in such systems there exists a pow- 
erful tool (Feshbach resonances) to control the strength 
of the interaction between atoms. Magnetically tuned 
scattering resonances allow for changing the magnitude 
as well as the sign of the scattering length which at 
low temperatures fully determines the atomic interac- 
tions. Recently, mixtures of fermionic gases turned out 
to be a successful way through in obtaining a molecular 
Bosc-Einstcin condensates (BEC) in thermal equilibrium 
. Such molecular condensates are good starting points 
to experimental study of BEC-BCS (Bardeen-Cooper- 
Schrieffer) crossover region; subject under intensive the- 
oretical investigation. Two-component Fermi gases en- 
abled also getting into a regime of strongly interacting 
fermionic systems Q, here confirmed by the observation 
of anisotropic expansion of the gas when released from a 
trap. In Ref. the pairing gap in a strongly interact- 
ing two-component gas of fermionic 6 Li atoms has been 
measured showing that the system was already brought 
in a superfluid state. 

Another idea of modifying the interactions between 
fermions is to immerse Fermi atoms into a Bose gas. In 
this way a degenerate gas of fermionic potassium ( 40 Ka) 
was forced to collapse by attractive interaction with ru- 
bidium ( 87 Rb) Bose-Einstein condensate |4j. As it was 
shown in this experiment, a sufficiently large number of 
atoms is required to bring a Bose-Fermi mixture to col- 
lapse. This result suggests that increasing the attraction 
between the bosonic and fermionic atoms one can induce 
the effective attractive interaction between fermions. In- 
creasing the effective fermion-fermion attraction could, 
on the other hand, lead to a superfluidity in fermionic 
component. This kind of superfluidity resembles what 
happens in superconductors where phonon-induced in- 
teraction between electrons becomes attractive. 

In Ref. 5] it has been shown (although in a one- 
dimensional case) that for strong enough attraction be- 
tween bosons and fermions a Bose-Fermi mixture enters 
a new phase where both Bose and Fermi components be- 
come effectively attractive gases. One of the character- 
istics of this regime is that it becomes possible to gener- 



ate bright solitons which are two-component single-peak 
structures with larger number of bosons than fermions. 
Since in a three-dimensional space such a mixture might 
show an instability that could destroy it, it is necessary 
to investigate in detail the stability-instability crossover. 
The structure and the instability of the boson-fermion 
mixtures have been theoretically studied in Ref. [6j, al- 
though the analysis was restricted to spherically symmet- 
ric systems with equal numbers of atoms in both compo- 
nents. The authors determine, by fitting the numerical 
data, the expressions for critical numbers of atoms as a 
function of scattering lengths. In this Letter we investi- 
gate the border of the stability region for any trap, find 
analytic formula for the critical radial trap frequency and 
compare it with results of numerical solution of hydro- 
dynamic equations describing the Bose-Fermi mixture. 

In the mean-field approximation the Bose-Fermi mix- 
ture of Nb bosons and Np fermions can be described in 
terms of atomic orbitals 0,0 that fulfill the following set 
of equations (j — 1, 2, Np) 
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and determine the many-body wave function of the mix- 
ture, which is of the form 
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We assume the contact interaction between atoms with 
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<7b and qbf being the coupling constants for the boson- 
boson and boson-fermion interactions, respectively. 

Eqs. Q imply that fermions move under the influence 
of the effective potential which is the sum of the harmonic 
trap and the potential that originates in the presence of 
bosons 
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In Fig. \I\we plot this potential (dashed line) as well as 
the fermionic density (the sum of orbital densities) in the 
case of one-dimensional system. It is clear from Q that 
when the bosons and fermions attract each other strongly 
enough, the second term in generates the well at the 
bottom of the harmonic trap (see Fig. 0. Hence, one 
can distinguish between fermions captured by the bosons 
(those with energies below zero) and others with higher 
energies that built broad fermionic basis. 
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FIG. 1: Fermionic density (solid line) normalized to one for a 
one-dimensional Bose-Fermi mixture of 1000 bosons and 100 
fermions. The strength of effective repulsion for bosons is 
qb = 0.0163 and the attraction between bosons and fermions 
is determined by qbf = —0.5 in oscillatory units. The dashed 
line shows the effective potential for fermions. 

To estimate the number of fermions pulled inside the 
bosonic cloud we use the Thomas-Fermi approximation. 
Since we consider the 87 Rb and 40 K atoms in a magnetic 
trap in their doubly spin polarized states, we utilize the 
relation mp(u;^ ) 2 = rafl(wf) 2 (with a equal to x,y, and 
z) and rewrite the potential J3J in the following form 
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where fig is the chemical potential for the Bose-Einstein 
condensate in the Thomas-Fermi limit 
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Here, a; lo = yJh/{m B ijJho) is the harmonic oscillator 
length and ujf- lo — {uixUJyLUz) 1 / 3 is the geometric aver- 
age of the oscillator frequencies. The scattering length 



as and the interaction strength g B are related through 
g B = 4Trh 2 a B /m B . 

Now, the number of fermions captured by bosons is just 
the number of states existing in the potential well defined 
by the upper branch of Since this potential acts in- 
side the bosonic cloud, one has to count only states with 
energies below ^f~IJ-B with respect to the minimum of 
the well. The well itself is the ha rmonic potent ial with 
frequencies multiplied by factors \/\g B F\/ gB + 1 in com- 
parison with original ones for bosons. The total number 
of states in a harmonic potential with the energy less 
than e is given by @ 
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Therefore, the number of fermions pulled inside the 
bosons equals 
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and after inserting the chemical potential given by (jSJ) 
this number becomes 
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Having calculated the number of fermions captured by 
bosons one can separate off the bosonic component. The 
energy of the Bose component of a Bose-Fermi mixture 
in the mean-field approximation reads 
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where 9 B = N B ^ B ' and n™(r) is the fermionic density 
drown into the bosonic cloud. Simplifying the descrip- 
tion of the system we introduce the Gaussian variational 
ansatz for the condensate wave function as well as for the 
fermionic density 
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where the radial and the axial widths, w and v respec- 
tively, are assumed to be the same for both species. This 
assumption allows for elimination of the fermionic de- 
grees of freedom. Now, the energy of bosons is given 
by the expression (in oscillatory units based on the axial 
frequency) 
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where 
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the aspect ratio (3 = uj^_/uj z , and the reduced mass 
fi = rriBrnF I (tub + mp). The first term in the formula 
(fTTfl is the kinetic energy of the Bose cloud, the second 
one describes the harmonic trapping potential whereas 
the last term covers the interaction between atoms. Since 
the bosons and fermions attract each other, cl B f is neg- 
ative and when this attraction is strong enough, i.e. in 
the region of parameters close to the stability border the 
constant A becomes negative. Therefore, on the edge of 
the mixture stability the expression Ijllfl can be treated 
just as the energy of a Bose gas of effectively attractive 
atoms. 

We now find, for a given trap, the maximum strength 
of attraction between bosons which still allows for the ex- 
istence of the stable Bose-Einstein condensate. We look 
for the minimum of the energy The necessary con- 

dition for that turns to the set of following equations for 
the widths 
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Solving the condition (|18|l determines the critical trap 
frequencies. It turns out that both terms on the left hand 
side of (|18|) are much bigger than the term on the right 
hand side. Therefore, the critical axial trap frequency is 
given by 
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The formula 11201) can be also regarded as a way of 
determining the mutual scattering length a B F when one 
knows the critical number of bosons (and at the same 
time the number of drown in fermions is given by (JSJ)) 
for a particular set of trap parameters. In such a case 
one has to solve the following equation 

CujI' 5 ^ |a Bi H 4 = (D \a BF \ + 1) 3/2 , (21) 
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Calculating w 2 from the first equation and inserting it in 
the second one leads to the equation 
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Since A is negative, the Eq. Ijl5|l has a solution when the 
maximum of the right hand side of (|15l) as a function of 
the width v is bigger than the value of the left hand side 
of l|15|l at the point of this maximum. Then the critical 
value of the strength A fulfills the biquadratic equation 
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generalizes well known result for a spherically symmetric 
trap 0- 

Equating (|12|) and l|17l) one obtains the condition for 
the parameters of the maximally compact trap which still 
holds the Bose- Fermi mixture of a given number of atoms 
and the interaction strengths 
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FIG. 2: Critical axial frequency as a function of the aspect 
ratio for the mixture of 10 4 bosons and 10 4 fermions. Points 
and lines come from the numerical solution of a set of hy- 
drodynamic equations and from the analytical formula 12UI . 
respectively. Solid circles (dotted line), squares (dashed line), 
and diamonds (solid line) correspond to the mutual scatter- 
ing length aBF = — 21.7nm, — 17.6nm, and — 13.8nm, suc- 
cessively. 

In Figs. [21 an d El w e plot v c ^ = ujf / (2ir) as a func- 
tion of the aspect ratio, calculated based on the formula 
l|20p. Since there is a controversy over the value of the 
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scattering length clef P- ITol ITlj. we included three dif- 
ferent values of it. We show also points obtained from 
numerical integration of a set of equations that are a 
hydrodynamic version of Eqs. The hydrodynamic 

equations can be derived from a set of equations for re- 
duced density matrices (see Ref. after making a 
local equilibrium assumption for fermions (i.e., utiliz- 
ing the Thomas-Fermi approximation) and calculating 
the interaction between bosons and fermions within the 
mean- field approach. The basic "objects" in the hydro- 
dynamic approximation are then the condensate wave 
function and the fermionic density and velocity fields. 
Hence, the equation for the condensate wave function 
looks like that in a set of Eqs. Q but with the last 
term equal to gsFriFf^ (with tif being the fermionic 
density). On the other hand, the equations describing 
fermions are the usual hydrodynamic equations, i.e., the 
continuity equation and the Euler-type equation of mo- 
tion. 
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FIG. 3: The same as in Fig. [5]but for larger number of atoms, 
here N B = N F = 10 s . 

Figs. and |3l show good agreement between the for- 
mula Q2U[) and the numerical points for weaker boson- 
fcrmion attraction (i.e., less negative scattering length 
clef)- The growing discrepancy for stronger attraction 
means that the number of fermions pulled inside bosons 
is overestimated. Indeed, the formula © derived based 
on the Thomas-Fermi approximation for bosons, does not 
account for the effect of shrinking the bosonic cloud while 
capturing more and more fermions. 

In conclusion, we have analyzed the Bose-Fermi mix- 
ture in a regime of parameters, where both gases start 
to behave as a systems of effectively attractive atoms 
and the existence of the mixture is in danger due to the 
possible collapse. We found an analytical formula which 
determines the parameters of maximally compact trap 
allowing for the stable mixture. We compare this esti- 
mate with the results obtained by solving numerically 
the hydrodynamic equations describing the Bose-Fermi 
mixture. It turns out that the collapse is very sensitive 
to the value of the scattering length asF and therefore 
could help to determine this scattering length. 
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